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Abstract. A fully-connected qubit network is considered, where every qubit interacts 
with every other one. When the interactions between the qubits are homogeneous, the 
system is a special case of the finite Lipkin-Meshkov-Glick model. We propose a natural 
implementation of this model using superconducting qubits in state-of-the-art circuit 
QED. The ground state, the low-lying energy spectrum and the dynamical evolution 
are investigated. Wc find that, under realistic conditions, highly entangled states of 
Greenberger-Hornc-Zeilingcr and W types can be generated. We also comment on the 
infiucncc of disorder on the system and discuss the possibility of simulating complex 
quantum systems, such as Shcrrington-Kirkpatrick spin glasses, with superconducting 
qubit networks. 
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Introduction 

The development of quantum information science has provided us with a fresh 
perspective on condensed-matter physics. In the last few years, new tools from quantum 
information theory have been applied to several problems in many-body physics [1] . At 
the center of this novel approach lies the problem of entanglement, that is, how to 
quantify the genuinely quantum correlations [2] in a many-body system and what these 
correlations can tell us about the system itself. The entanglement in a given system 
is also a resource for various quantum information processing tasks. Spin chains and 
lattices with short-range interactions have been studied extensively in this context. 
Although interactions between neighbors are more common, there are instances where 
long-range interactions give a better description of physical systems, such as certain 
types of spin glasses [3], H]. It would be very interesting to be able to use quantum 
networks as simulators of such complex quantum systems. Furthermore, connected 
networks are attracting considerable interest in quantum information science [5] and in 
many other fields [6]. 

Here we study a fully-connected network, where every qubit interacts with every 
other one, irrespectively of the distances between them. Crucially, the proposed model 
is readily implementable with superconducting qubits [7] . The emerging field of circuit 
QED [71 El El [in] provides a natural system in which a large number of qubits can be 
coupled together. In such systems, superconducting qubits play the role of atoms, and a 
harmonic-oscillator circuit element plays the role of a cavity with which they interact. If 
a single 'cavity' is simultaneously coupled to a number of qubits, it will mediate coupling 
between all the possible qubit pairs [11] . If, in addition, the cavity is far off resonance 
with the qubits, its degrees of freedom can be integrated out of the problem and we 
obtain a system in which all the qubits are pairwise interacting. Previous studies have 
considered similar circuits for coupling arbitrarily distant superconducting qubits |12j . 
However, these studies relied on time-dependent pulses to selectively couple one qubit 
pair at a time, whereas here we consider the simultaneous coupling of all qubit pairs. 
An important incentive for studying the fully-connected network is that all the different 
elements for its construction are already in place in the laboratory. 

Next we introduce the model, study its low-lying energy spectrum and its dynamical 
response, and discuss the influence of disorder. Our analysis mainly concerns the 
entanglement properties of small networks and the generation of highly-entangled states 
in near-future experiments with existing technologies of superconducting qubits in 
circuit QED. We also discuss the possibility of simulating spin glasses with these systems. 

Model and Hamiltonian 

To begin with, we consider charge (flux) qubits that are coupled capacitively 
(inductively) and assume that each qubit is operated at its degeneracy point [7]. 
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Therefore the Hamiltonian is 

^ A, 

'H= -TT^i ~ Jij^i^j^ (1) 

where the second sum runs over all possible qubit pairs. Here, > is the level 
splitting and Jij is the strength of the coupling between qubits i and j. Xi, Yi and Zj 
denote the Pauli matrices for qubit i. Multi-qubit entanglement generation has been 
analyzed in trapped ions and atoms using Hamiltonians of the form in Eq. ([T]) [T3l [T^ . 
In a circuit QED setup, which is the focus of this work, the macroscopic qubits allow 
individual addressing and readout, in addition to relatively straightforward scalability. 
Although additional terms will appear in the Hamiltonian of this system, these terms 
can be made negligibly small under realistic conditions, as was shown in Ref. [TT] . 

If we let Aj = A and Jij = J for all qubits, the system is homogeneous and it 
corresponds to a special case of the Lipkin-Meshkov-Glick model [H, [15]. In this case, 
the Hamiltonian can be expressed as 

where .^Totai = Y^iLi^i and Xxotai = Sili^i- From Eqs. ([1]) and ([2]) it is clear 
that the Hamiltonian commutes with the square of the total pseudo-spin operator, 
J2i{Xi + Yi + Zi)^ , and possesses spin-flip symmetry, i.e., it also commutes with HjZj [15]. 



Ground state properties 

Two parameter regimes can be identified, namely, A ^ A^|J| and A <^ A^|J|. In the 
first case, the single-qubit term in Ti dominates and the preferred basis is the eigenbasis 
of Z, {|0), |1)}, with Z\0) = |0) (the eigenstates of X are denoted by {|+), |-)}). As 
I J| ^ the ground state of the system becomes equal to the fully-separable state 

I^SEP) = |0)^^. (3) 

We shall employ a pseudo-spin language for convenience; when Jij < we say that the 
interaction is 'antiferromagnetic' (AFM), and when Jij > we say that the interaction 
is 'ferromagnetic' (FM). 

Secondly, if A ^ N\J\ then the interaction term in Ti, dominates and the preferred 
basis is the |— )} basis. In this case the sign of J becomes important. For large, 
positive J (FM regime) the interaction term tends to set all the qubits in the same state 
in the {|+), |— )} basis. The 'ideal' ground state of the system is approximately 




which is known as the GHZ-state |2]. In practice, however, this ideal state is typically 
fragile under small external perturbations. By contrast, for large and negative J (AFM 
regime) the interaction term favors a state in which pairs of neighbouring qubits are 
antiparallel in the {|+), |— )} basis. Clearly, in a fully-connected geometry this condition 
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Figure 1. Fidelity |(^'g|0)| between the ground state, |^'g), and the states I^'sep), 
I^GHz), I^'ent) against J for a network of iV = 3 qubits. We show (a) the ideal case 
and (b) the case where a small perturbation gXi of strength g = A/ 100 is applied to 
qubit i. 



is impossible to satisfy because every qubit neighbors every other qubit. This high degree 
of frustration in the system leads to highly-entangled states, as we show below. For the 
particular case of = 3 there is a relatively simple ground state, namely, 

I^ent) = ^(| + -0) + |-0+) + |0 + -)). (5) 

We illustrate the above statements by means of an exact numerical diagonalization, 
Ti.l'^n) = -E'nl^n) (for n = 1, 2, . . .). In all numerical simulations, we let A = 1 and hence 
express the results in units of A. In Fig. [l](a) we show the fidelity |(\E'g|0)| between 
the ground state of the system, I^E^g); and the three states 10) of Eqs. ([3]) to ([5]) for 
= 3, in the ideal case. In Fig. [11(b) we calculate the same fidelities but now we apply 
a small perturbation gXi, of strength (7 ^ A, to qubit i (it does not matter which 
one). The symmetry-breaking term only affects the FM regime, where the I^E'ghz) 
becomes increasingly fragile as J is increased beyond a certain optimal value. This 
result holds for all small networks with ~ 10: it is always possible to find an optimal 
value of the coupling strength such that the ground state is very close to a GHZ-state, 
even in the presence of a small external perturbation. For larger networks, or stronger 
perturbations, the behaviour is more abrupt and we do not obtain exact or almost exact 
GHZ-states, in practice. 

In Fig. [2](a) the full energy spectrum for the N = 3 case is presented. In Fig. [2](b) 
the lower half of the spectrum for = 4 is presented. One can see that the ground 
state is unique in the AFM regime and two-fold degenerate in the FM regime. Also, 
there is a finite energy gap between the ground states and the first excited states, which 
remains constant with increasing \J\ in the AFM regime and increases with \J\ in the 
FM regime. 

We now turn to the entanglement properties of the different possible ground states. 
We consider both the entanglement between pairs of qubits and also between blocks of 
qubits. Among the various measures of two-qubit entanglement [2], we calculate the 
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logarithmic negativity [16]. It is defined as 

^N(pij) = iog2 IIpJII, (6) 

where 1 1 . 1 1 denotes the trace norm of a matrix and pjj is the partial transpose of the 
reduced density matrix pij of two qubits i and j. We also use the von Neumann entropy 
of a state pk, for a block of k < N qubits, 

'S'(pfe) = -tr(pfclog2Pfe)- (7) 

S{pk) quantifies how mixed is the reduced density matrix pk and, if the system as a 
whole is in a pure state, it also quantifies the amount of entanglement between the 
qubits in the set k and those in the rest of the system. 

The pairwise entanglement, E^{pij), is shown in Fig. [3](a) against J for systems of 
size = 3, 4, . . . , 12. We observe that E^{pij) decreases with the size of the system TV. 
This can be explained by the fact that as increases so do the number of interactions 
for each individual qubit. This higher degree of connectivity places stronger monogamy 
constraints that generally reduce the two-party entanglements [17] . To study the many- 
body correlations of the ground state, we show in Fig. [3t^b) the entanglement of a single 
qubit i with the rest of the system, using the von Neumann entropy S{pi). In the FM 
regime the ideal ground state is a GHZ state, in which every qubit is maximally entangled 
with the rest of the network. In practice, however, under the influence of an external 
perturbation, the system chooses one of the two degenerate eigenstates and becomes 
fully separable in the deep FM regime. This is seen in the inset of Fig. [SlJ^b), which 
shows the ground state of a small network perturbed by gXi with strength g = A/100. 
In the AFM regime, S{pi) increases with the size of the network for odd and it 
decreases for even A^. This is due to the different symmetries in the two cases, which 
appear for small networks but disappear in the thermodynamic limit (see, e.g., [15j). 
In summary, in the AFM regime the system achieves multi-qubit entanglement for any 
I J| > 0; in the FM regime the system achieves a high degree of multi-qubit entanglement 
in practice, approaching a GHZ-state, for an optimal value of the interaction strength, 
while deep in the FM regime it becomes fully separable. 
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Figure 3. (a) E-[<i{pij) against J for networks of size iV = 3 to = 12 (from top 
to bottom), (b) S{pi) against J for = 3 to = 10. Solid lines correspond to 
A^ ~ 3,5,7,9 (from bottom to top) and broken lines to A^ = 4,6,8,10 (from top to 
bottom), as indicated. In both (a) and (b), the insets correspond to the N = 5 case 
with a small perturbation gXi of strength g = A/100. 



Dynamical evolution 

Next, we study the system's evolution, which is determined by the state 

|^(t)) =exp(-zm)|^(0)). (8) 

We consider two simple initial states: the state |\&sep) of Eq. ([3]), and the state 
I^SEp) = |1) ® lO)*^^""^, which is the same as I^E'sep) except that the state of one 
qubit is flipped. Both of these states are separable and easy to prepare experimentally. 
For definiteness we assume that J > 0, and we only consider the case of weak coupling, 
which is most relevant to near-future experiments. 

In the weak-coupling limit, J ^ A/A^, one can classify the low-lying states 
according to the number of elementary excitations they contain. The ground 
state I^E'sEp) contains no excitations and does not evolve in time. At an energy 
~ A above the ground state, there are energy eigenstates that can be 
identified as one-excitation states. These states have the 'spin- wave' form |$i,a:) = 
Ef=i e^'^^^'^/^lOi ■ ■ ■ Oj_iljOj+i ■■■On), where k = 0, 1, ■ ■ ■ , - 1. Their energies 
(relative to the ground state) are given by A — (A^ — ^)J for k = and A-\- J otherwise. 
As these A^ energy eigenstates are separated from all other states by an energy at least 
~ A, their dynamics can be analyzed in the restricted Hilbert space containing only these 
A^ states. Using the above spectrum, we find that the initial state |\E'(t = 0)) = |^sep) 
evolves, up to an overall phase factor, into 

I^W) H^I'^P) + ^"P''^" ^ v .). (9) 



where the 'W-state' is \Wn) = Ef=i |0i ■ ■ ■ Oj_iljOj+i ■■■On)- From this result we 
observe that the evolved state is a time-dependent superposition of I^'sep) IVTat). 
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For large networks (A^ — * oo) the second term can be neglected, and the system remains 
close to its initial state (the initial excitation is localized). For = 3 and = 4, on 
the other hand, Eq. ([9]) reduces to variants of the W-state (when |e*^"^* — 1| = ^/N), 

„iti7r/6 p±5i7r/6 

m = ^|100) + -^yf- (1010) + 1001)) 

IW4) = i(-|1000) + lOlOO) + lOOlO) + lOOOl)). (10) 

For larger, but still weak, coupling we study the evolution using the fidelity between 
|\Ef(t)) and the GHZ and W states. If the initial state is |\1/sep), then the system evolves 
into a state that is close to a GHZ-state; if the initial state is I^I^sep)' then the system 
evolves into a state that is close to a W-state, for small networks. For = 3 and 
N = 4, (t)) can be made arbitrarily close to |\&ghz) or |VV) for properly chosen values 
of J (see Fig. Hj). It is possible to achieve multi-qubit entanglement of these two types 
with networks of size A^ > 4, but the window for the required coupling strengths is 
much sharper and the fidelity maxima decrease below 1. For instance, for A^ = 6 and 
< J < A/2 the maximum fidelity of the evolved state with the GHZ-state is 0.96, 
which is still a relatively high fidelity [T8] . 

Thus, the fully connected network is well suited for the fast (one-step) preparation 
of GHZ- and W- entangled states. This result should be contrasted with recent work 
on the generation of entangled states in superconducting qubit circuits using generally 
long sequences of basic operations [19] . For networks with A^ > 4 and the simple initial 
states used above, we obtain entangled states that are different from the GHZ- and 
W-states, but still highly entangled and therefore of potential value for future quantum 
information technologies (e.g., we only mention applications related to quantum secret 
sharing and quantum average estimation [20]). 

Finally, we note that the dynamics of entanglement has simple periodicity only 
in the case of A^ = 3, which corresponds to a closed chain. In general, for A^ > 3, 
the behavior of the logarithmic negativity E^lpij) is complicated and there are time 
intervals for which it drops to zero. As the system size increases, so does the durations 
of these time intervals. During most of the intervals in which the pairwise entanglements 
vanish, the block entanglements increase to their local maxima. Therefore the system 
achieves a form of multi-qubit entanglement. On the other hand, for small and negative 
J we find that the evolution of entanglement shows more periodic features. The block 
entropies S{pk) (for different block sizes k) have their local maxima at the same times 
as the EnIpij). Hence in the AFM regime the system evolves into an entangled state 
that is different from the GHZ- and W-states. 

Static disorder 

In the presence of disorder, a system with the Hamiltonian 7Y(Aj, Jij) of Eq. (fl]) is 
inhomogeneous. As a result, different partitions of the network that correspond to the 
same number of qubits are no longer equivalent. We study the influence of uniform, 
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Figure 4. Fidelity between the evolved state and the different initial states, 

the GHZ-state of Eq. (g]) and the W-state of Eq. ((10|), as indicated, (a) and (b) 
correspond to = 3, J = A/2, while (c) and (d) correspond to = 4, J = A/3. 



static disorder on the system by assuming that the Jij and the Aj are chosen randomly 
from the intervals [1 — 6i,l + 6i]J and [1 — ^2, 1 + 62] A, respectively, where 61^2 G [0, 1] 
quantify the amount of disorder in each parameter. We neglect the type of disorder 
that can move a qubit away from its degeneracy point, but this is a valid approximation 
for realistic implementations. The properties of the low-lying energy sector and the 
dynamics are studied as an average over many realizations of ?i(Aj, Jj^). We only 
consider the weak-coupling regime, and the case of small networks. 

We shall not show any details on these Monte-Carlo simulations, as they do not add 
any insight beyond the main conclusions. From the numerical calculations, we observe 
that disorder in the Jij is more important than disorder in the Aj. More crucially, we 
find that the results reported here, including those on the ground state entanglement 
and the dynamics of entanglement, remain largely unaffected for disorder of amount 
<^i,2 < 10%. This result is in agreement with previous studies on disorder [21], and it is 
an experimentally achievable upper bound. 

Spin glasses 

The simplest examples of glassy systems are spin glasses and they offer the possibility 
of studying the behavior of complex systems away from equilibrium. Spin glasses 
arise when the interactions between spins are ferromagnetic for some bonds and 
antiferromagnetic for others, in which case the spin orientation cannot be uniform in 
space even at low temperatures [3] . In this case the spin orientation can become random 
and frozen in time. A particularly illuminating and extensively studied model of spin 
glasses is the Sherrington-Kirkpatrick (SK) model. The SK model in a transverse field 
[1] is given by the Hamiltonian Ti of Eq. ([T]) for interactions that are disordered in both 
magnitude and sign. In fact, if the distribution of Jij is Gaussian, then for J ~ A the 
system is a spin glass for temperatures lower than ksT ~ A/4 (see, e.g., [1]). 

There are various important open problems in this field, both theoretical and 
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experimental. For instance, in relation to the proposal of the present work, one such 
problem is the behavior of spin glasses in the low temperature phase region, where 
quantum phenomena dominate (see, e.g., Parisi in [3]). Some of the advantages of using 
qubit networks as quantum simulators include the fact that the states of all the qubits 
can be prepared controUably and that the dynamics of all the qubits can be monitored 
as a function of time, yielding the 'microscopic' dynamics of individual spins and not 
just averaged spin quantities. Therefore the implementation of fully-connected networks 
with superconducting qubits, which can be addressed and measured individually, can 
offer valuable additional tools for the study of complex quantum systems, such as spin 
glasses. 

Outlook and summary 

In this work we have focused on XX-type interqubit couplings [Eq. ([T])] since this is 
relevant to present-day experiments in circuit QED. A modified version of the flux 
qubit was proposed recently, which implements ZZ-tjpe couplings p2]. This form of 
the coupling could give rise to other types of many-body entangled quantum states |23j . 
We have also focused on the case where the cavity is far off resonance with the qubits. 
Bringing the cavity into resonance with the qubits would lead to a star geometry with 
the cavity at the center of the star. In this case, the Hamiltonian of Eq. (fl]) is no 
longer valid; instead, each qubit interacts only with the cavity and not other qubits. 
The addition of a nonlinearity to the cavity can be used to make the center of the star 
an effective two- level system [2^ . 

In conclusion, we have proposed an implementation of a qubit network where all 
qubits are coupled in pairs, independently of the relative distances between them, as 
in the finite LMG model of spin systems. Such a network supports a highly-entangled 
ground state in the case of antiferromagnetic interactions, which is robust against small 
external perturbations. Under suitable conditions, separable initial states evolve into 
exact GHZ- and W-type states in the case of small networks, or other highly-entangled 
states for larger networks. The qualitative behavior of the system is unaffected by 
the presence of static disorder, as long as the amount of disorder is under about 10%. 
Thus the system is well-suited for the generation of many-body entangled states with 
macroscopic superconducting qubits. The presence of entanglement in such systems can 
be witnessed experimentally via combinations of two-qubit correlation measurements, 
as for example described in Ref. [21]. Another promising prospect is the simulation of 
spin glasses. The fully-connected network could be realized experimentally in the near 
future with superconducting qubits in circuit QED. 
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